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Abstract 

We make a quenched lattice simulation of hadron spectroscopy at j3 = 6.2 
with the Wilson action non-perturbatively improved. With respect to the unim- 
proved case, the estimate of the lattice spacing is less influenced by the choice 
of input hadron masses. We study also the effects of using an improved quark 
mass in the fits to the dependence of hadron masses upon quark masses. 
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1 Introduction 



The computational cost of the extrapolation to the continuum limit of lattice QCD 
simulations can be significantly reduced by using improved actions, where the leading 
cutoff effects are cancelled by suitable counterterms. It has been shown that 0(a) 
on shell improvement is achieved by adding to the usual Wilson action the clover 
term, with a coefficient that has been recently determined non-perturbatively by the 
ALPHA collaboration @, §. 

We present here the results of our study of hadron spectroscopy using the non- 
perturbatively clover-improved Wilson action 

S = S w + c sw a 5 - £*r» a„ v F^(x) , (I) 

X 

where Sw is the standard Wilson action and the second term is the naive continuum 
limit of the clover term (see @ for details). Preliminary results have been reported 
in 0. 

We consider a lattice volume of 24 3 x 48 and a coupling (3 = 6.2. According to 
0, [|, we thus take csw = 1.61375065. We choose the following values for the hopping 
parameter k: 0.1240, 0.1275, 0.1310, 0.1340, 0.1345, 0.1350, 0.1352. The simulations 
were carried out on the 512 processor computer of the APE100 series at Tor Vergata. 

Our statistics come from 104 quenched gauge configurations, generated by a hy- 
brid over-relaxation algorithm, with each update corresponding to a heat-bath sweep 
followed by three over-relaxation sweeps. The configurations are separated by 1000 
updates. 

The inversion of the fermion matrix is performed using the stabilized biconjugate 
gradient algorithm Q. We restart the inversion from the current solution every 100 
iterations, in order to reduce the accumulation of roundoff errors ||. We employ 
point-like sources. We sum fermion propagators over the space directions x, y, z for 
sites within blocks of side 3, and then store the result. We then form hadron cor- 
relations from these "packed" propagators. This procedure differs from the exact 
computation of hadron correlations by gauge-non-invariant terms and becomes exact 
in the limit of an infinite number of configurations (in our case we have checked that 
the errors thus introduced are negligible). This corresponds to gaining a factor 3 3 in 
storage, and has enabled us to have all quark propagators stored simultaneously. In 
this way we can form "off-line" hadron correlations from non-degenerate combinations 
of quark flavours. 

Hadron masses are obtained from single mass fits to the large time behaviour 
of zero momentum hadron correlators. A two mass fit has also been done, but the 
results for the higher mass were too unstable to quote numbers. The value for the 
lower mass turned out to be totally compatible with the single mass fit. The errors are 
estimated through a jack-knife procedure. In order to improve the stability in time of 
the plateau where a single hadron dominates the correlation for the baryons, in some 
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cases we have averaged the correlation over a fixed number of lattice spacings in the 
time direction: the resulting correlation maintains its exponential behaviour with a 
coefficient depending upon the size of the smearing. This procedure turned out to be 
very useful for the determination of baryon mass splittings where we also made fits 
directly to the ratio of correlations to minimize the effects of collective fluctuations. 
We report in Table [I] our results for hadron masses in lattice units and at the various 
values of k. 



K 


M PS 


My 


M N 


0.1240 


1.0737(20) 


1.0997(15) 


1.7002(35) 


0.1275 


0.8532(20) 


0.8873(15) 


1.3774(35) 


0.1310 


0.6048(20) 


0.6565(20) 


1.0189(40) 


0.1340 


0.3445(20) 


0.4364(35) 


0.6608(65) 


0.1345 


0.2909(25) 


0.3977(50) 


0.5937(85) 


0.1350 


0.2294(30) 


0.364(10) 


0.520(15) 


0.1352 


0.2007(40) 


0.353(15) 


0.483(20) 



Table 1: Masses in lattice units (diagonal-flavour combinations only) for the pseu- 
doscalar and vector mesons, and for the nucleon. 



2 Lattice Spacing and Meson Masses 

In order to extract physical values from lattice data, one needs to interpolate the 
results at appropriate values of quark masses, and to give them physical values by a 
suitable normalization. These steps imply a number of choices. 

The first concerns the determination of "k c ": indeed, the renormalization of the 
quark mass with Wilson fermions is not multiplicative, and the critical value of the 
hopping parameter is shifted from its free-case limit. Within the improvement pro- 
gramme one can determine the critical value k c using the mass extracted from an 
improved Ward identity, defined as 

(d^A^ + c A ad^ bare )}0) 
mwi = 2(p( bare )0) ' 

with the parameter ca fixed from ref. to —0.037. A linear extrapolation to the 
limit of my/ 1 = , using the first five points obtained from the combinations of the 
three highest k values, provides a fit for the determination of n c much more stable 
than the conventional fit of pseudoscalar meson masses to the limit of zero pion mass. 
Our results for K c £1X6! 

from my/! = k c = 0.135828(5) 
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from M 2 p S = k c = 0.135849(13). 
The values of the parameters of a linear fit in 1/k to Mp s and to m-^j are: 

m WI = -3.8059(15) + 0.5169(2) x (1/k) (2) 
M PS = -8.38(3) + 1.138(4) x (1/k). (3) 

The residual discrepancy between the two values of k c is a signal of residual lattice 
artefacts. 

Once the choice for the value for k c from mwi has been made, we make a second 
choice that improves the fits to the dependence of hadron masses upon quark masses: 
instead of using the bare quark mass, defined as m q (K) = (1/k — 1/k c )/2, we use an 
"improved" bare quark mass @ defined by 

m g (K) = m,(«) [1 + b m m q (K)} . (4) 

(Note that m q is the renormalized mass with Z m = 1 .) The improvement coefficient 
b m has been determined non-perturbatively @: b m = —0.62(3) . 

For non-degenerate flavour cases we use symmetric averages of the masses de- 
fined above, e.g. for a meson corresponding to the flavours k± and k 2 , we define 
m g (Ki,K 2 ) = [m 9 (/ci) + m ? («2)]/2. Similarly, for baryons, we define m q and m q as 
the symmetric average of the quark masses for the three flavours. Note that when 
using averages of the improved masses we reabsorb in the mass definition terms that 
are quadratic in the quark masses. This may and does change the estimates of the 
deviation of hadron masses (their squares for pseudoscalars) from a linear dependence 
upon the quark masses. For baryons the choice of the improved mass does extend the 
linear behavior substantially beyond the light-mass region. 

A third choice to be made is on the physical inputs that give physical units to the 
lattice spacing and to the lattice masses. 

The lightest of our quark masses is in the strange quark mass region and we decided 
to use physical inputs at the strange quark mass for the lattice spacing to avoid the 
inclusion of systematic uncertainties deriving from a chiral extrapolation. A known 
deficiency of the quenched approximation is the estimate of vector-pseudoscalar mass 
splittings. In particular the experimental difference of the squares of vector and pseu- 
doscalar masses is essentially constant as a function of the quark mass, a feature that 
is not reproduced by present quenched data. Besides possible unquenching effects, 
this can be ascribed to lattice artefacts, which should be partially cancelled in our 
simulation. 

For the strange mesons, we decided to use as input just such a square- mass dif- 
ference: this can be seen as an attempt to reabsorb the residual lattice artefacts in 
a redefinition of the quark mass. More in detail, we first fix the value of the strange 
quark mass from the ratio: 

(M 2 -M 2 PS )/M 2 (5) 
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using the experimental value for K and K* meson masses as input, and then we 
extract the lattice spacing by normalizing the K* mass to its experimental value. 
This gives the following values: 

m s = 0.0315(45) and a" 1 = 2561(100) MeV , (6) 

where we have used, for the average values, quadratic and linear mass fits, respectively. 
These are the values of m s and a -1 adopted in this paper. In the following we will 
not include the overall error coming from the uncertainty of the lattice spacing when 
quoting errors for hadron masses in physical units. This amounts to an overall and 
systematic uncertainty of about 4%. 

Our fits for My and Mp s in the light/strange quark mass region are given below. 



• Linear fits in terms of the improved quark mass: 

M v = 0.309(20) + 2.49(55) x m q (7) 

M 2 PS = -0.001(2) + 2.398(35) x m q (8) 

• Quadratic fits in terms of the improved quark mass: 

My = 0.298(15) + 2.82(50) x m q + 1.0(35) x m\ (9) 

M 2 PS = 0.005(2) + 1.992(50) x m q + 7.53(40) x m\ (10) 

• Linear fits in terms of the unimproved quark mass: 

M v = 0.311(20) + 2.38(55) x m q (11) 

M% s = 0.000(2) + 2.328(30) x m q (12) 

• Quadratic fits in terms of the unimproved quark mass: 

M v = 0.304(15) + 2.60(45) x m q + 0.4(30) x m q 2 (13) 

M 2 PS = 0.004(2) + 2.085(50) x m q + 3.91(35) x m 2 (14) 



The quadratic fit for the pseudoscalar meson gives evidence for a quadratic term, 
while for the vector case its presence appears very marginal. The coefficient of the 
quadratic term in the unimproved mass case is expected to be smaller than in the 
improved case from the relation in eq. @. 

The chiral limit of vector masses is obtained by linear extrapolation in terms of the 
improved quark mass, showing good agreement when compared to the experimental 
value. This indicates that a lattice spacing extracted from the chiral limit leads to a 
value very close to the one in eq. @: 

a: h \ ral = 2486(175) MeV. (15) 
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Figure 1: APE plot for My. 

The fit for the vector case goes through the <fi mass with reasonable accuracy and 
indicates that in the strange quark mass region a lattice spacing normalization from 
the (f) would also give similar results: 

a^ 1 = 2634(40) MeV. (16) 

In Fig. |l] we show a plot of the vector meson mass as a function of Mp s . The 
masses are normalized to the K* meson mass, taken at the strange quark mass value 
given in the eq. (^j above. The asterisks correspond to the physical values for the 
ratio. 

As can be seen from this figure, the experimental slope dMy / dMp s is different 
from the one corresponding to our data. This is related to the behaviour of the 
quantity J defined by: 

J = My {dMy/dM 2 ps ) , (17) 

obtained from My as a function of Mp s . Our results for J are quoted in Table |[ 
They are rather similar to those obtained without improvement (for a recent review, 
see H), although slightly closer to the experimental result. 

The procedure used to determine the strange quark mass fails when applied to the 
charm case: it is not possible here to find a value of the quark mass that reproduces 
the observed splitting. Therefore, we decided to define the charm quark mass from a 
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Exp. value 


Our value 


Jk* 


0.487 


0.40(11) 


J(f> 


0.557 


0.45(12) 



Table 2: Values of the quantity J and comparison with experiment. 



fit to the ratio M 2 PS jM 2 K ,. We get 

m c = 0.388(25) from m q , (18) 
m c = 0.3154(10) from m q . (19) 

By using the U K*" lattice spacing we obtain for the D* charmed meson: 

M D * = 1955(6) MeV (exp. value 2008 MeV). (20) 

The above value for the vector mass comes from a linear fit to a cluster of points 
around the improved charm quark mass. 

For the vector meson with an unimproved quark mass we find: 

M D * = 1981(5) MeV. (21) 

In both cases, the underestimate of the experimental pseudoscalar-vector mass split- 
ting is not dramatic. 

The value of the lattice-improved strange quark mass given in eq. @ can be 
used to obtain the strange quark mass in physical units. Firstly, we multiply it by 
the mass renormalization factor Z m at the (1/a) scale, using the tadpole- improved 
formula given in eq. (41) of ref. ||. We then evolve the result to the 2 GeV scale with 
a perturbative renormalization group factor as in [IU|. We get 

m s = 111(15) MeV . (22) 

This result is in agreement with that of ref. ||. 

In the ratio of charm over strange quark mass the renormalization factor Z m drops 
out and we can give a value for the ratio of renormalized quark masses: 

m c/ m s = 12.1 ± 1.6 from m q , 
m c /m s = 10.0 ± 1.4 from m q . 

Notice that the difference comes essentially from the reduction of the charm quark 
mass going from the unimproved to the improved bare quark mass definition (see eq. 

We can use this ratio, combined with the theoretical prediction of 1.220(60) GeV 
for the charm quark mass |TTJ, to obtain an independent evaluation of the strange 
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quark mass. After evolving this value of the charm quark mass to the 2 GeV scale, 
we get: 

m s = 92(13) MeV from m q , 
m s = 111(16) MeV from m q . 

The value obtained using the improved bare quark mass is in very good agreement 
with the determination in eq. (p2|) . 

3 Baryon Masses 

In this section we report our results for baryon masses and baryon mass splittings. 
We show APE plots for the octet and for the decuplet baryons in Fig. |2|. Note that we 
divide at each point by My, interpolated to the strange-quark mass, corresponding 
to Mk* (i.e. a constant value). Experimental points in these figures (asterisks) cor- 
respond to M?v, My, and Mn and appropriate meson masses in Fig. 0a, and similarly 
M A , .Uv-..\/h and M n for Fig. |b. 

We have considered non-relativistic wave functions. The correlators for the spin- 
3/2 (decuplet) baryons are completely symmetric in flavour. We therefore expect to 
observe a smooth behaviour of their masses in terms of the average quark masses m q 
and m q . For the spin-1/2 (octet) baryons, on the contrary, the correlators are not 
completely flavour-symmetric. We have two types of correlators (see [BJ for details): 
S-like (e.g. for the proton, neutron, E and S) and A-like (e.g. for A). The expressions 
for their masses from quenched chiral perturbation theory, in the case of two light 
quarks m u and a strange quark m s , are given by |J: 

M s = M + 4Fm u + 2(F - D)m s (23) 
M A = M + 4(F-2D/3)m tt + 2(F + D/3)m s . (24) 

The constant D can be related to the S-A octet hyperfine splitting, as can be seen 
from the above formulae. We have considered fits as functions of the symmetric 
average masses described above also in the spin-1/2 case. 
Note that we have defined the quantity 

M oct = (Mv + M A )/2 , (25) 

which we construct from the appropriate ratios of correlators. This combination was 
chosen so that the resulting mass is flavour-symmetric. For this "octet" particle the 
advantage of using the improved quark mass is clearest, see Fig. [| 

From fits of the E-A mass splitting, and of M oct , it is possible to extract the value 
of the constants D and F. We obtain: 

D = -0.51(20) (26) 
F = 0.80(20). (27) 
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Figure 2: APE plot for (a) the nucleoli mass (octet baryons) and (b) the A mass 
(decuplet baryons). 
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Figure 3: "Octet" mass in terms of (a) m q (b) m q . 
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In Table [3] we give our baryon mass values in MeV. We have used in all cases linear 
fits of the masses as functions of the average improved bare quark mass m q , except 
for the E-A mass splitting, where we used a quadratic fit. We have also included the 
predictions for charmed baryons using the improved quark mass. 





Exp. value 


Our value 


M N 


939 


952(110) 


M E 


1189.4 


1148(60) 


M E _ A 


73.7 


70(30) 


M A 


1232 


1265(110) 




293 


297(80) 


M Sc 


2453 


2406(15) 




168 


232(40) 



Table 3: Baryon masses in MeV and comparison with experiment. The error in 
parentheses represents the statistical error, while the error due to the determination 
of the inverse lattice spacing is of about 4%. 



4 Conclusions 

The main effect of the non-perturbative improvement observed on hadron masses 
is a smaller spread in the value of the lattice spacing extracted from light mesons 
with and without strange quarks. We have also shown how the use of the improved 
quark mass turns the rough behaviour of the dependence of the octet baryon mass 
upon quark masses into a smooth one. We have found a remarkable agreement in 
the improved theory between two independent determinations of the strange quark 
mass, one normalized through the lattice spacing and the other from the value for 
the charm mass extracted in the continuum from charmonium spectrum calculations. 

For an accurate comparison of these results with the unimproved case, we are 
presently analysing ]12|] on the same set of gauge configurations the case c sw = 0. 
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